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Theory of nonlinear cyclotron resonance in quasi-two-dimensional electron systems 
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Momentum and energy balance equations are developed for steady-state electron transport and 
optical absorption under the influence of a dc electric field, an intense ac electric field of terahertz 
(THz) frequency in a two-dimensional (2D) semiconductor in the presence of a strong magnetic 
field perpendicular to the 2D plane. These equations are applied to study the intensity-dependent 
cyclotron resonance (CR) in far-infrared transmission and THz-radiation-induced photoconductivity 
of GaAs heterostructures in Faraday geometry. We find that the CR peaks and line shapes of the 
(in ' transmittance exhibit different intensity dependence when the intensity of THz field increases in 

^ , the range above or below a certain critical value. The CR in photoresistivity, however, always 

enhances with increasing the intensity of the THz field. These results qualitatively agree with 
' the experimental observations. We have clarified that the CR in photoconductivity is not only 

, the result of the electron heating, but also comes from photon- assisted scattering enhancement, 

especially at high temperatures. The effects of an intense THz field on Faraday angle and ellipticity 
of magnetically-biased 2D semiconductors have also been demonstrated. 

o . 

^ , PACS numbers: 73.50.Jt, 73.50.Mx, 78.67.Dc, 78.20.Ls 

I. INTRODUCTION 

' Cyclotron resonance (CR) is a fundamental process of carriers in quasi-two-dimensional semiconductors subjected 
simultaneously to a magnetic field and a far-infrared or terahertz (THz) ac field. It occurs when the separation 
I , between two adjacent Landau levels is closed to the photon energy of the ac field, and leads to the resonant behavior 
in absorption and in photoconductivity as functions of the magnetic field. 

CR in absorption or transmission has been proved to be a powerful tool providing valuable information of two- 
^ ' dimensional (2D) semiconductors on the carrier scattering processes, the confining potential, and even the spin 
relaxation mechanismsi Generally, weak far-infrared irradiations were used in most of experiments, and theoretical 
analyses focused on the linear response of the system to the ac fieldr^ When the strength of the incident radiation 
field increased, nonlinear behavior of CR in absorption and transmission was observed. It was found"^ that the 
. resonant field to zero field transmissivity ratio descended slightly with increasing the intensity of the far-infrared field 

■ from 0.1 W/cm^, then increased for radiation-field intensities above lOW/cm^. This kind of intensity-dependent CR 
lO transmission has not yet been explained so far except a simple comparison with linear-response formula of Drude-type 
OO taking the scattering time as a fitting parameter. 

. Another manifestation of CR shows up in photoresponse, namely, photoconductivity or photoresistivity. The latter 

■ is defined as the longitudinal dc magnetoresistivity change induced by the irradiation of the high-frequency electro- 
■»«^ ' magnetic wave, which exhibits a resonant peak structure under CR condition. This effect has long been known at 

low temperatures, and was believed to arise from the electron heating induced redistribution of the electron gas upon 
absorbing the radiation-field energy^i^ The energy absorption and thus the rise of the electron temperature exhibit 
'J^ maxima when the photon energy of the radiation field equals the cyclotron energy, yielding the resonance peak struc- 
'■^ ture of the photoresistivity. Later analysis of experimental data^ indicated that, in addition to this heating-induced 
^ , electron redistribution, another nonthermal mechanism might be needed to account for the observed photoconductiv- 
O ' ity of GaAs/AlGaAs heterostructures at CR. Recent photoresistivity measurement of GaAs/AlGaAs heterojunction 
performing at lattice temperature T — 150 K subjected to irradiations of 4 THz frequency,^ surprisingly showed re- 
markable CR peaks. At such a high lattice temperature with strong polar optic phonon scattering providing an 
efficient energy dissipation channel, the radiation-induced electron-temperature increase is far smaller to account for 
such strong CR in photoresistivity. So far no convincing nonthermal mechanism responsible for far-infrared photo- 
conductivity in semiconductors has been proposed. A nonlinear theory for high-temperature photoresistivity CR in 
two-dimensional system is still lacking. It is an urgent need to develop a tractable, microscopic theory capable of 
treating nonlinear absorption and photoresponse under intense THz radiation and to clarify the missing nonthermal 
mechanism for high-temperature photoresistivity in two-dimensional semiconductors. 

A few years ago, one of the authors developed a balance-equation approach^ for hot-electron transport driven by 
a THz electric field of single frequency (E^sm{ujt)). This method made use of the fact that, when the harmonic 
generation is small and the frequency gets into the THz regime or higher, the electron drift velocity in the steady 
transport state oscillates almost out of phase of the electric field, i.e. the drift velocity is essentially of the form 
Vi cos(wi). At the same time, all orders of this (frequency w) photon assisted impurity and phonon scatterings are 
included in the relaxation processes. This method has been successfully applied to discuss THz photoabsorption and 
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THz-induced dc conductivity response in bulk and two-dimensional semiconductors in the case without magnetic 
fieldSiS or with a magnetic field in Voigt configuration^ However, this assumption of such time-dependent form of 
drift velocity will no longer be true when there is a strong magnetic field B not parallel to Vi. Since the Lorentz 
force evi x Bcos([i;<) acts on the moving electron, its steady-state drift velocity will contain a term {ijJc/lo)vi sin(a;t), 
where Uc — \eB\/m is the cyclotron frequency. This velocity, which is perpendicular to vi and B, and oscillates 7r/2 
out of phase of vi cos(cL;i), is of the same order of magnitude as vi in the frequency and magnetic field range u) ~ ujc- 
Because of this, the balance-equation method used in Ref.Hol is not able to deal with the problem related to cyclotron 
resonance in semiconductors in Faraday geometry. 

The purpose of this article is to develop a theory for electron transport driven by an intense far-infrared electric field 
in a quasi-2D semiconductor in the presence of an arbitrary magnetic field perpendicular to the 2D plane. We find that 
it is possible to extend balance equation method proposed in Ref.|3 to the case in the presence of a quantized magnetic 
field in Faraday configuration, and to deal with those problems in magnetotransport, such as the intensity-dependent 
THz transmission and photoresponse at CR in quasi-2D electron systems. 

This paper is organized as follows. In Sec. II we will sketch the derivation of the force- and energy-balance equations 
of quasi-two-dimensional electron systems subjected to a magnetic field, an arbitrary dc and THz ac electric fields. In 
Sec. HI the cyclotron resonance in time-dependent drift velocity is discussed. The investigation on cyclotron resonance 
in transmission and photoconductivity is, respectively, presented in Sec. IV and Sec. V. Finally, in Section VI a short 
conclusion will be given. 

II. FORMULATION 
A. Hamiltonian 

We consider iVe electrons in a unit area of a quasi-two-dimensional system, such as a hetero junction or a quantum 
well. In these systems the electrons are free to move in the x-y plane, but are subjected to a confining potential V{z) 
in the z-direction. These electrons are interacting with each other and also coupled with phonons and scattered by 
randomly distributed impurities in the lattice. 

To include possible elliptically polarized electromagnetic radiation wc assume that a uniform dc electric field Eg 
and a terahertz ac field E(t) of angular frequency tj, 

E(i) = Es sin(cjt) + E^ cos(wt), (1) 

are applied in the x-y plane, together with a uniform magnetic field B = (0, 0, B) along the z axis. These magnetic 
and electric fields can respectively be described by a vector potential A(r) and a scalar potential (p(r, t) of the form 

V X A(r) = B, (2) 

^(r,i) = -r-Eo -r •E(t). (3) 
In the presence of these electric and magnetic fields the Hamiltonian of the system has the form 

H = HeEit) + Hei + Hep + Hph ■ (4) 

Here 

+ '^^cirq -rjn,z^,zj) (5) 

i<j 

is the Hamiltonian of electrons driven by the electric and magnetic fields with Vc being the electron-electron Coulomb 
interaction, Hd and -ffop are, respectively, the electron-impurity and electron-phonon couplings. In equation |(SJ, 
r^ii = {xj,yj) and Pj|| = (pjxiPjy) are the coordinate and momentum of the jth. electron in the 2D plane, and z and 
Pjz are those perpendicular to the plane; m and are, respectively, the effective mass parallel and perpendicular to 
the plane. 

It is convenient to introduce two-dimensional center-of-mass momentum and coordinate variables P = {P-x, Py) and 

R = {Rx, Ry)- 



p = Ep.ih 




1 

2m 
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and the relative-electron momentum and coordinate variables = (p'^pPjz) and r' = (r^n , Zj){j = 1, ...Nc): 

1 _ 



R. 



(7) 



In term of these variables, the Hamiltonian iJ^E can be separated into a center-of-mass part -ffcm and a relative 
electron part Her- 



Hp, 



Hem ^ , 



2Nem 

E 



(P - iVceA(R))2 - TVoeEo • R - A^ceE(t) • R, 



' 'p;ii-^A(r;,|);' 



2m 



(8) 
(9) 

(10) 



It should be noted that the relative-electron Hamiltonian Ha is the just that of a quasi-2D system subjected to a 
uniform magnetic field in the z direction. Its eigenstate can be designated by a subband index s, a Landau level index 
n and a wavevector ky, having energy spectrum 



£sn = £s + - ^)wc, s = 0, 1,... and n = l,2,. 



(11) 



where Uc = \eB\/m is the cyclotron frequency. In this paper we will limit to the case that the 2D electrons occupying 
only the lowest subband and ignore index s. In Eq. (5) the electron-impurity and electron-phonon interaction H^i 
and i?ep have the same expression as those given in Ref. llll in terms of center-of-mass coordinate R and the density 
operator of the relative-electrons, 



(12) 



B. Balance equations 

On the basis of the Heisenberg equation of motion we can derive the velocity (operator) of the center-of-mass, V, 
which is the rate of change of the center-of-mass coordinate R, and equation for the rate of change of the center-of-mass 
velocity V = R: 

V--i[R,i/] = -^(P-7VeeA(R)), (13) 
iVem 



and 



with 



-i[V,i/] + ?X = ^{Eo+E(0+VxB} + -^, (14) 



F = -i^[/(q||,z,)q||e''iii-(^-"^")pqi, -i^Af(q,A)0qAe"ir%q||. (15) 
We can also derive equation for the rate of change of the relative electron energy: 

Tier = -i[i?er,i?] = "i ^ C/(q|| , 0a)e''»" (^"'■"^Pq,, - i ^ A/(q, A)0qAe''l" '^Pq,, . (16) 
q||,a q,A 

In Eqs. (|15I16|) {Ya,Za) and C^(q||, Za) are the impurity position and its potential, Af (q, A) is the matrix element due to 
coupling between electrons and a phonon of wave vector q = (q|| ,qz) in branch A having energy fiqA, (/"qA = 6qA + ^-qA 
stands for the phonon field operator, and pqn = — i[pq|| , iJcr]- 

In order to derive the force- and energy-balance equations we need to carry out the statistical average of these 
operator equations. As proposed in Ref. 12, we treat the center-of-mass coordinate R and velocity V classically, and 
by neglecting their small fluctuations we will regard them as time-dependent expectation values of the center-of-mass 



4 



coordinate and velocity, R(i) and V(t). In the present paper, we are mainly concerned with the steady transport 
state under a terahertz irradiation of single frequency and focus on the photoresponse of the dc conductance and the 
energy absorption and transmission of the THz signal. These quantities are directly related to the time-averaged and 
base-frequency oscillating components of the center-of-mass velocity. The second and higher harmonic components of 
the electron velocity, if any, though give no direct contribution to ac-field transmission and dc photoresponse, would 
enter the frictional force, energy transfer and energy absorption rates in the resulting equations, thus may in turn 
affect the time-averaged and lower harmonic terms of the drift velocity. However, unless for a specially designed 
system, the generated power of the third harmonic current in an ordinary semiconductor is generally less than a few 
percent of the base-frequency power even in the case of strongly nonlinear transport when the ac field amplitude gets 
as high as several kV/cm or higher Its effect on the photon-assisted scattering matrix element will be an order 
of magnitude further weaker than that of the base frequency photonsiSi24i& For the radiation field intensity concerned 
in the present study, which is an order of magnitude smaller than that for the above-mentioned harmonic generation, 
the effect of higher harmonic current is safely negligible. Hence, it suffices to assume that the center-of-mass velocity, 
i.e. the electron drift velocity, consists of only a dc part vq and a stationary time-dependent part v(i) of the form 

= Vq + vi cos(wt) + V2 sm{ujt) . (17) 

On the other hand, in order to carry out the statistical average we need the density matrix p. For 2D systems having 
electron sheet density of order of, or higher than, 10^^ m~^, the intrasubband and intersubband Coulomb interaction 
are sufficiently strong that it is adequate to describe the relative-electron transport state using a single electron 
temperature T^. Hence, the density matrix can be obtained from solving the Liouville equation by starting from the 
initial state of the relative-electron-phonon system at time t = — oo, in which the phonon system is in equilibrium at 
lattice temperature T and the relative-electron system is in equilibrium at an electron temperature Tc'. 

p\t=-oo = Po = le-^"/^<=e-^p^/^ (18) 

Z is the normalized parameter. 

With the density matrix thus obtained to the first order in Hd + Hep, we can carry out the statistical average of 
operator equations H14() and 116|1 . In the procedure the form factor 



rt °° 



VQ — nuj){t — t') 



71— — OQ 



m/0 



Q — nuj){t — t') 



n— — oo 



(19) 



is expanded in terms of the Bessel functions J,i(a;). Here ^ = -\/ (q|| ■ vi)^ + (q|| . m^)^ juj and satisfies the relation 
tan(^ = (q • V2)/(q • vi). 

Since all the transport properties are measured over a time interval much longer than the period of terahertz field and 
we are concerned only with photoresponse and photoabsorption of the system, it suffices for us to know the frictional- 
force for the time oscillating term with base frequency lo and the energy-related quantity for the time-averaged term. 
The frictional force can be written as 

F(t) = Fo - (Fii - F22) sin(u;t) + (F12 + F21) cos(tji), (20) 

with the functions Fq and F^^, (/x, v = 1,2) given by 

00 oc 

Fo=^q|||{/(q||)|' 4(e)n2(q||,^o-riL^) + 5Iq||l^'^(q'^)l' E J'(OA2(q,A,coo + r!qA-na;), (21) 

qy 71— — 00 q,A n— — 00 



Fi^ = -5]q||?7p|t^(q||)|' E K(0]'ni(q||,a;o-nu;)-5]q||r7p|M(q,A)|' ' Ai(q, A, c^o + 17^^ -nc^), 

qy n— — 00 q, A n— — oo 

(22) 



F2^=Eq||^|(7(q||)|' 2n4(e)H2(q||,c.o-n^)+EqilTl^^('l'^)l' E 2nj2 (e)A2(q, A, c^o + f^q^ 



71 — — 00 



71 — — OC 



(23) 
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In these expressions, 77^ = qy - Vf^/uj^; luq = qy -Vq; f/(q||) is related to the impurity potential U{q\\,Za) and distribution 
of impurities along z axis A2(q, A, fl) is the imaginary part of the electron-phonon correlation function, which 

can be expressed through the imaginary part of electron density correlation function 112 (q|| , ^) 



A2(q,A,n) = 2n2(q||,fi) 



(24) 



with n(x) = l/[exp(x) — 1] being the Bosc function. The real parts of the electron-phonon correlation function and the 
electron density correlation function, Ai(q, A, fl) and ni(q|| , f2), can be obtained from their imaginary parts through 
Kramers-Kronig transformation. 

The momentum-balance equation obtained by taking the statistical average of operator equation H14|l . has the 
following form 

1 e 

ViLo smiuit) - V2UJ cos(ujt) = — — F(t) + — {Eq -I- E(i) -I- [vq + v(t)] x B} . (25) 

m 

That is 

= A^eeEo + iVee(vo x B) + Fq, (26) 
VI = — - (Fn - F22) - — (V2 X B), (27) 

eE 1 e 

-V2 = — + - (F12+F21) (vixB). (28) 

muj I\cmuj mio 

The energy-balance equation is obtained by taking the long-time average of statistically averaged operator equation 
(HHJ to be 

TVecEo -^0 + Sp-W = 0. (29) 

Here W is the time-averaged rate of the energy transfer from the electron system to the phonon system, whose 
expression can be obtained from the second term on the right side of equation (|21l) by replacing the qy factor with 
i^qA- 5'p is the time- averaged rate of the electron energy gain from the radiation field and has the following form 

00 00 
5p = ^|C/(q||)|' riu;Jl{On2{qii,u;o-ncj) + Y,\M{<l,X)\^ ^ ncc-J^ (^)A2(q,A, c^o + f^qA - "cc>). (30) 

q|| n— — C30 q.A n— — oo 

Note that Sp is negatively equal to time averaged Joule heat {v{t) ■ F{t))^ = iVee(Eo • Vq -I- Eg • V2/2 -I- Ec • Vi/2). 

Momentum and energy balance equations H26() to (|29|l constitute a close set of equations to determine the parameters 
Vq, Vi, V2, and Te when Eq, Ej, and E^ are given. 

The sum over n in the expressions for Eg, F^^, W and Sp represents contribution of all orders of multiphoton 
processes related to the photons of frequency w. In the present formulation, the role of the single- frequency radiation 
field is two fold. (1) It induces photon-assisted impurity and phonon scatterings associated with single (|n| — 1) and 
multiple (|ri| > 1) photon processes, which are superposed on the direct impurity and phonon scattering {n = 0) term. 
(2) It transfers energy to the electron system (Sp) through single and multiple photon-assisted process. 

Note that n2(q||,f2) and ni(q||,i7) are respectively the imaginary and real parts of the electron density correlation 
function of the 2D system in the presence of the magnetic field. In the Landau representation, one can write^ 

n2(q||,r!) = -i2^C„,„.(^''?y/2)n2(n,n',r!), (31) 



71, n 



n2{n,n',n) = / de[/(e)-/(e + f7)]ImG;(£ + n)ImG;;,(e), (32) 



where I = ^l/\eB\ is the magnetic length. 



Cn,n+i{Y) = ^y'e-^[LUr)]2 (33) 
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with Ll^{Y) being the associate Laguerre polynomial, f{e) = {exp[(£ — /x) /Te\ + 1}~^ is the Fermi distribution function, 
and ImGJ^(e) is the imaginary part of the Green's function of the Landau level n, which is proportional to the density 
of states, such that the density of electrons is given by 

= E / rf£/(£)IniG;(£). (34) 

This equation determines the chemical potential. 

In principle, to obtain the Green's function of Landau levels n, G^{t), a self-consistent calculation has to be 
carried out from the Dyson equation for the self-energy with all the scattering mechanisms included. ^8, The resultant 
Green's function is generally a complicated function of the magnetic field, temperature, and Landau-level index n, also 
dependent on the relative strength of the impurity and phonon scattering. In the present study we do not attempt 
a self-consistent calculation of GJ^(i). Instead, we choose a Gaussian- type function for the Landau- level shape for 
simplicity^ 

G;;(0 = -ze{t) exp[-*(n - ^)u;J - ^Tlf] (35) 

with a unified broadening parameter r„ = F for all the Landau levels, which is taken as [2euJc/i^'rn^QY/'^ . When the 
hot-electron effect is neglected, /io corresponds to the linear mobility at temperature T in the absence of magnetic 
fieldsi^- In order to consider the hot-electron-induced Landau levels broadening, we will empirically treat /ig as the 
linear mobility of the system in the absence of the magnetic field at temperature T^. Note that this Gaussian- type of 
Green's function is proved to be correct at low temperature^ On the other hand, it also has been used to interpret 
the magnetophonon resonance at high lattice temperature^ In the present paper we will show that this Green's 
function can lead to a qualitative agreement between theoretical and experimental results within the magnetic field 
range considered. To improve the agreement further, a more careful study on the Green's function of a Landau level 
should be performed. 

Above formulation can be used to describe the transport and optical properties of magnetically-biased quasi-2D 
semiconductors subjected to a dc field and a terahertz field. The conventional magneto-optical study in the far- 
infrared frequency regime corresponds to the case of zero dc field Eq = 0, where one studies the intensity-dependent 
terahertz absorption, transmission and other effects in the presence of a strong magnetic field. On the other hand, to 
investigate photoconductivity, we should treat the weak dc field limit of our formulation. 

III. CYCLOTRON RESONANCE IN DRIFT VELOCITY 

Substituting the force Eq. (|27|l into Eq. (|28(l , we can write 



e 



e 



1 /_ _ e 



vi = (1 - w.Vw^)-^ E, + (E,xB) Fn-F22 [(F12 -t- F21) x B] U , (36) 



Nf,moj V moj 



V2 = (l-c.2/a,2)-i| — 



Ee (E, X B) 



— ^ (Fi2 + F21 + — [(Fn - F22) X B] 



(37) 



Cyclotron resonance is easily seen in the case of weak scatterings when the terms with ¥ functions in the above 
equation are small: both vi and V2 exhibit peaks at CR. Since all the transport quantities, including W , Sp and Fq, 
are functions of the drift velocity as well as the electron temperature To, and the latter is determined by the energy 
balance equation (30), the CR of vi and V2 will result in the CR in W, Sp, Fq and To. 

Eqs. H36|l and (|37|1 can be further simplified when the radiation field is weak and the dc field is absent. In this case 
can be treated as small parameters. To the first order of these small parameters, the force function F^^ can be 
written as 



(38) 



where Af^(w,vo) are the real (/i = 1) and imaginary (/i — 2) parts of the memory functionsii* It is convenient to 
write out the expression for the complex velocity u+ = vix + iw2x and u_ = V2y — Wiy rather than for Vi and V2 



eT/2m* 
—er /2m 



-B4 



{ijJ—uI)t + \ (cLl+CJ*)T + i 

E— 



{lO—UJ*)t + 1 {LJ + LO*)T + i 



(39) 
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Here, we have defined 

m* = to[1 + vo)/a;], 

1/r = Af2(co,vo)/[l + Mi(cj,vo)H, 

Lo* ~ eB/jn*, 

E+ ~ Esx + Ecy + i{Esy — Ecx), 

E- — Esx — Ecy — i{Esy + Ecx), 

with (Esx^Esy) = 'Es and (Ecx, Ecy) = E^. 

Our weak-field results for vi and V2 reduces to those of Ref._2 in the case of circularly polarized ac fields. 



IV. TRANSMISSION 



For normally incident electromagnetic wave, the transmitted electric field E(i), which is regarded as the field driving 
the 2D electronSfSi is related to the incident electric field Ei(i) by 

, Ncev(t)/eoc 2nn „ , , 
E(t) = — — ^ " H ^Ei(t). (40) 

Us + no ns + no 

Here no and Us are the relative refractive indices of the airs and 2D semiconductors, and c and eo are the light speed 
and the dielectric constant in vacuum, respectively. The transmitted field E(i) depends on the drift velocity v(t) of 
the 2D system. In the following numerical studies on transmission and photoconductivity, we will assume a sinusoidal 
incident field Ei(i) = (E'i^ sin(a;t), 0) along x axis and derive E(i) self-consistently together with v(t). 

We have numerically calculated the magneto-optical properties of a GaAs/AlGaAs hetero junction subjected to 
a THz ac field and a magnetic field. The strength of transmitted ac field E(t), and the parameters v^(/x=l,2) 
and Tc, are obtained by resolving the Eqs. H4()|l . H27|l . (|28|l and (|29|l . We consider a GaAs-based quasi-2D system 
having electron density = 2.5 x 10^^ ni^^ and 4.2 K linear mobility 50m^/Vs (which is used to determine the 
impurity density) at lattice temperature T =4.2 K, similar to that used in Ref. 3. The elastic scattering due to 
randomly distributed charged impurity and the inelastic scattering due to polar optical phonons (via Frohlich coupling 
with electrons), longitudinal acoustic phonons (via deformation potential and piezoelectric coupling), and transverse 
acoustic phonons (via piezoelectric coupling with electrons) are taken into account. The material and electron-phonon 
coupling parameters are taken as typical values for GaAs. In the numerical calculation the maximum Landau level 
is taken to be 20, and the summation over multiphoton indices n are carried up to a given accuracy of 10^^ for each 
quantity. 

The transmittance T, defined asSi 

T^iml2L (41) 

with < .. >t denoting the time average. When connecting it with measured quantities, of course, the multiple 
interference between interfaces of the substrate has to be taken into acount. The calculated transmittance and 
the corresponding electron temperature are plotted in Fig. 1 as functions of the intensity of the THz field at two 
frequencies uj/2tt = 0.83 and 1.6 THz in the center position of CR, namely ojc = It can be seen that, the 
transmittance first decreases gently with increasing intensity of the THz radiation from zero, and reaches a bottom at 
a critical intensity around lOW/cm^, then increases rapidly with further increasing the field strength. This feature 
appears more pronounced at lower frequency, in consistence with the experimental observation,"^ as shown in the inset 
of Fig. 1, where the measured transmittance for 1.6 THz and 0.24 THz exhibits similar trend. The case of 0.83 THz 
does not show minimum. This deviation is believed to come from the experimental errors. In Fig. 2 we display the 
transmittance CR line shape for incident electromagnetic fields of different intensities at frequency 0.83 THz. The line 
width exhibits no significant change below the critical intensity but increases rapidly when the intensity of THz field 
grows above the critical value. 

This kind of i?is-dependent behavior of transmittance is in agreement with the intensity dependence of the absorption 
rate a ~ Su/Ef^ in the absence of the magnetic field. For a similar 2D GaAs-based semiconductor without magnetic 
field, Ref.|^ showed that, the absorption percentage increases with increasing strength of the radiation field from 
low-field value, then reaches a maximum (of order of 2 percent) at the field amplitude around several kV/cm before 
decreasing with further increase of the radiation field strength, and that lower frequency has stronger maximum. 
At low velocity side, when hot-electron effect is relatively weak and the direct impurity and phonon scatterings 
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change little, the behavior of the absorption rate comes mainly from the drift- velocity dependence of the multi-photon 
assisted scattering matrix element, as described by the Bessel functions J^(^) in the expressions for Sp. In fact, all 
the multi-photon {n > 1) contributions to the absorption coefficient are zero at vanishing velocity and reach maxima 
at finite (increasing with n) drift velocities, the resultant absorption coefficient first increases with increasing velocity. 
When the drift velocity becomes sufficiently large, reduction of absorption rates, induced by the large argument of 
the lowest-order Bessel functions, will exceed the increased contributions from the other multi-photon processes. This 
leads to the drop of the absorption rate. In the present case having a strong magnetic field, CR greatly enhances the 
drift velocity vi and V2 at Wc ~ for a given Eis in comparison with the case without magnetic field. Therefore, the 
maximum absorption rate should appear at much smaller strength of radiation field and have much larger value for 
the case of cyclotron resonance than in the absence of a magnetic field or far away from CR. 

The behavior of the transmittance CR line shape is related to Landau level broadening due to hot-electron effect. 
For THz field below the critical intensity the electron temperature is less than 60 K (see Fig. 1) and impurities are 
the dominant scatterers, yielding almost a constant mobility fj,o (thus the Landau level broadening). When the THz 
field goes above the critical intensity, the electron temperature grows rapidly and polar optical phonons become the 
dominant scatterers, giving rise to a strongly temperature-dependent mobility /io, thus a Landau level broadening 
which increases rapidly with increasing field strength. 

Our formulation can also be employed to investigate the change of the transmitted electromagnetic field polarization 
in quasi-two-dimensional electron systems. This effect is known as Faraday effect and has long been investigated under 
linear condition. The present approach provides an convenient formulation to study the Faraday effect for the 
case when the incident light is strong and the nonlinear absorption occurs. 

The relevant quantities characterizing the Faraday effects are the ellipticity ry and Faraday rotation angle 0p , which 
are determined through the amplitudes of transmitted field E(t) 

tan?/ = {a+ - a^)/{a+ + a^), (42) 
ep = (0+-0-)/2. (43) 



Where 



tau0+ = |f^L±|£^, (44) 
tanr = 'f-+^-\ (45) 



= J{E,.,+E,yY + {E,,~EsyY/2, (46) 



_ [Es, - E,yy + {E,, + E,yf/2, (47) 

with [Esx.Esy) = Es and {E^x^Ecy) = E^. 

We plot the calculated results of -q and 6-p in Fig. 3 for the above mentioned two-dimensional sample. The resonance 
in ellipticity and antiresonance in Faraday angle can be seen evidently. Their line shapes also manifest different 
behavior when the intensity of THz field lies below or above the critical value. These intensity-dependent behaviors 
of ellipticity and Faraday rotation can also be understood by multi-photon-assisted scatterings and hot-electron effect 
induced Landau level broadening. 

V. PHOTOCONDUCTIVITY 

The response of the linear dc conductance to far-infrared irradiation is easily obtained in the the weak dc field limit 
of our formulation. Taking vq to be in the x direction, Vq = (uo£!;,0,0) and expanding the equation (|26|) to the first 
order in vqx, we obtain the transverse and longitudinal resistivities Rxy and Rxx as follows: 

Rxy = = B/N,e, (48) 



^n.(q,„i^) 



i1— — nuj 



^ oo 



^A.(q,A,l^) 



(49) 

O— OqA —nuj 
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The parameters Vi, V2 and To in these expressions should be determined by solving equations (|27() . (|28l) and H29|) 
with zero Vq. The longitudinal photoresistivity is defined as 



(50) 



with being the longitudinal magnetoresistivity in the absence of the radiation field. 

Photoconductivity in semiconductors, in the absence or in the presence of magnetic fields, has long been known 
at low temperatures, and was understood to result from the effects of electron heating due to the absorption of the 
radiation field energy^*SS*^ In our formulation, the photoconductivity arises not only from the hot-electron effect 
(electron temperature change), but also from the photon-assisted electron-impurity and electron-phonon scatterings. 
Although it is difficult to distinguish contributions to photoconductivity from different mechanisms when the applied 
terahertz field is strong, in the case of weak ac fields, the longitudinal photoresistivity can be written as the sum of 
two terms: 



(51) 



The first term Ai?£^ is obtained through expanding Eq. (|49|l by the small parameter ATc = Tc — T and is the result 
of ac field induced electron temperature change, as that proposed first by Kogan for the case without a magnetic 
field^ After determining the small electron temperature change from energy-balance equation (|29|l , we can write 



Here 



and 



n 



AR^j!} = $ • AT,. 



2^g^|M(q,A)| — ^n( — ) 



8 5]^|Af(q,A)pn2(q||,l]q,)n'(%^ 



|yn.(q,„i7) 



(52) 



(53) 



y2 



T 



(54) 



with 



N,mLoM2iio,0) ~ J2 f^qAqji |M(q,A)|'A2(q||,nqA±c^) 

q,A,± 



(55) 



Note that the terms on the right hand side of Eq. H55|l are respectively the changes of Sp and W induced by photon- 
assisted scatterings. 

The second component of photoresistivity AR. 
field induced electron distribution change 



?(op) 
XX 



Ai?l°P^ 



1 



4Af2e2, 



is the result of the photon-assisted scattering processes and ac 
(3<+3«L + <+0Q2(^). (56) 



Where 



Q2(^) ^J2lt [^<i(q|hO) -^rf(q||,c^)] /iJ, 



(57) 



with 



A,(q||, fl) - |C/(q||)|' ^n2(q||, f]) + J] |M(q, A)|' J^U^, A, n). 



g.,A 



(58) 



In the weak ac field limit we can see from Eq. H39() that the amplitudes of time-dependent drift velocity are linearly 
dependent on the strength of applied THz fields. Consequently, the photoresistivity is proportional to the intensity of 
THz field. When the intensity of the driving field becomes strong, the dependence of photoconductivity on the strength 
of THz field exhibits a complicated behavior. At the same time, the contribution from the electron-temperature 
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change and nonthermal photon-assisted scattering are hybridized. At low lattice temperatures, the hot-electron effect 
is sufficiently strong and is generally the dominant mechanism for photoresistivity. At high temperature, however, 
when the polar optical phonon scattering provides an efficient energy dissipation channel, the photoconductivity is 
mainly contributed from nonthermal mechanism. 

Recently, the photoconductivity CR at high temperature has been demonstrated in the experiment of Ref. 0- 
The remarkable peaks in photoconductivity are observed when the cyclotron frequency is closed to the frequency of 
terahertz fields. It is also found that, not only the height but also the width of CR peaks increases with increasing 
intensity of terahertz fields. 

In order to illustrate CR in photoconductivity at high lattice temperatures, we have numerically evaluated the dc 
longitudinal photoresistivity of a magnetically-biased GaAs/AlGaAs heterojunction. The lattice temperature T is 
150 K. The considered sample has electron density ~ 2.0 x 10^^ m~^ and 4.2 K linear mobility 200m^/Vs, similar 
to that used in the experiment of Ref. 0- 

In Fig. 4 the longitudinal photoresistivity induced by THz fields of frequency uj/27t — 4 THz having several different 
amplitudes is plotted as a function of magnetic field. The resonant structure near the cyclotron resonance position 
shows up clearly. Furthermore, with increasing strength of THz field, the CR peaks ascend and the line shapes 
broaden. These features are in qualitative agreement with experimental results in Ref.Q In the inset of Fig. 4 we plot 
the photoresistivity at CR as a function of the intensity of THz field. One can see that the photoresistivity follows 
a linear dependence on the intensity of THz field in the range < lis < IkW/cm^. For larger ac field intensity, the 
deviation from linear dependence appears. The electron temperature also exhibits resonant peak when the cyclotron 
frequency is closed to the THz frequency, as shown in Fig. 5. Nevertheless, since at high lattice temperatures the 
strong electron-LO phonon scattering provides an efficient energy dissipation channel, the rise of electron temperature 
is modest, and the hot-electron effect induced photoconductivity is small in comparison with the nonthermal effect. 
The photon-assisted scatterings are the main mechanisms for the photoconductivity at high temperature. For relative 
weak radiation field, e.g. Eis < 0.24kV/cm, the electron temperature has no appreciable difference from TjVet 
the ARxx still exhibits sizable resonance, comparable with the peak height observed in experiments of Ref. In 
Fig. 6 we plot the longitudinal photoresistivity as a function of magnetic field strength for relatively weak ac fields. 
The contribution to photoconductivity from hot-electron effect is less than 3 %. We also show the experimental 
observation^ of CR in photoconductivity in the inset of Fig. 6. The extraordinary width of the resonance shown in the 
experimental AR^x is possible due to the short pulse or large bandwidth of the terahertz radiation used in experiment, 
as explained in Ref. 

VI. CONCLUSION 

We have developed the momentum and energy balance equations for steady-state electron transport and optical 
absorption under the influence of a dc electric field, an intense THz ac electric field in a two-dimensional semiconductor 
in the presence of a strong magnetic field perpendicular to the 2D plane. This formulation allows us to investigate the 
THz-field-intensity dependence of the cyclotron resonance in transmittance and photoconductivity of GaAs/AlGaAs 
hetero junctions. We found that the CR peaks and line shapes of transmittance exhibit different behaviors when the 
intensity of the THz field increases in the range above or below a certain critical value. The cyclotron resonance in 
photoresistivity, however, always enhances with increasing intensity of the THz field. These results qualitatively agree 
with the experimental observations. The intensity-dependent behavior of transmittance at CR is explained as the 
result of combing hot-electron effect induced Landau level broadening and electron-phonon scattering enhancement, 
and the drift velocity dependence of the photon-assisted scattering matrix elements. We have also clarified that the 
CR in photoconductivity is not only the result of the electron heating, but also comes from photon-assisted scattering 
enhancement, especially at high temperatures. The effect of an intense THz field on Faraday angle and ellipticity of 
magnetically-biased 2D semiconductor systems have also been demonstrated. 
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Figure Captions 

FIG. 1. The intensity-dependence of 2D semiconductor transmittance and electron temperature, is plotted at 
cyclotron resonance position, lUc — uj. The transmittance is normalized to the value of zero magnetic field. The THz 
fields with two different frequency uj/2Tr = 0.83 and w/27r = 1.6THz are exposed to the studied 2D system. The 
lattice temperature is T = 4.2 K. Experimental results of transmittance versus the intensity of THz fields lis (Fig. 2 
of Ref. is reproduced in the inset. 

FIG. 2. The cyclotron resonance of 2D semiconductor transmittance of several incident electromagnetic fields with 
a same frequency uj/2tt = 0.83 THz but different intensity lis ^ 0.1, 0.6, 74, 750, 1500, 2500 W/cm^. 

FIG. 3. The Faraday angle 9^ and ellipticity rj are plotted as functions of the strength of magnetic fields B for the 
same system and under the same condition as described in Fig. 2. 

FIG. 4. The cyclotron resonance in the longitudinal resistivity change AR^x induced by a radiation field of frequency 
cj/27r ~ 4THz having several different amplitudes Eis — 0.24, 0.61, 0.87, 1.1, 1.2 kV/cm. The lattice temperature is 
T — 150 K. The inset shows AR^x versus intensity of incident THz fields lis aX ujc — ^■ 

FIG. 5. The cyclotron resonance in electron temperature To for the same system and under the same condition as 
described in Fig. 4. 

FIG. 6. Similar to the Fig. 4, but the strength of THz fields are relatively weak Eis = 0, 0.086, 0.13, 0.19, 0.24kV/cm. 
The experimental results (Fig. 3 of Ref. 0) of the cyclotron resonance in photoconductivity is shown in the inset. 
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